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We investigate the interaction induced superfiuid to Mott insulator transition in the one- 
dimensional (ID) Bose-Hubbard model (BHM) for fillings n = 1, n = 2, and n = 3 by studying 
the single-particle gap, the fidelity susceptibility, and the amplitude of Bloch oscillations via density 
matrix renormalization group (DMRG) methods. We apply a generic scaling procedure for the gap, 
which allows us to determine the critical points with very high accuracy. We also study how the 
fidelity susceptibility behaves across the phase transition. Furthermore, we show that in the BHM, 
and in a system of spinless fermions, the amplitude of Bloch oscillations after a tilt of the lattice 
vanishes at the critical points. This indicates that Bloch oscillations can serve as a tool to detect 
the transition point in ongoing experiments with ultracold gases. 
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Introduction — Ultracold atomic and molecular gases 
in optical lattices provide a unique playground for inves- 
tigating quantum many-body phenomena [1, 2]. Since 
the seminal experiment by Greiner et al. [3], it has be- 
come common in such experiments to study quantum 
phase transitions in the presence of strong correlations. 
In particular, optical lattice realizations of the BHM 

H--Ji: {blb^ + H.c.) + I E". - l) ' (1) 
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have been shown to undergo a transition from a super- 
fluid to a Mott-insulator as the ratio of C// J is increased 
in different dimensions [3-5]. In what follows, we set 
J = 1 and ft = 1, so that U is measured in units of J and 
time t is measured in units of h/ J. 

The ID BHM, the focus of this study, is of particular 
interest because of the dominant role played by quantum 
fluctuations. From the theoretical side, it is challenging 
to accurately determine the critical value Uc at which the 
system at constant density undergoes a superfluid-Mott- 
insulator transition, something that, due to the lack of 
exact solutions for this model, is typically done utiliz- 
ing computational approaches [2]. Here, the Berezinskii- 
Kostcrlitz-Thoulcss (BKT) universality class of the tran- 
sition makes calculations in finite systems susceptible to 
large finite-size effects. Insights from Luttinger liquid 
theory, combined with DMRG [6-9] calculations of cor- 
relation functions and extrapolations to the thermody- 
namic limit, have provided some of the most accurate 
values of Uc to date [10, 11] (see Ref. [2] for a review). 
Due to the large numerical effort needed, alternative and 
more accurate scaling approaches to calculate Uc, which 
do not rely on computing correlation functions, are highly 
desirable. From the experimental point of view, many of 



the quantities used in theoretical studies to determine Uc 
are either difficult (e.g, the gap [4, 12-14]) or not possi- 
ble to measure accurately. The task is complicated even 
further by inhomogeneities induced by the unavoidable 
confining potentials present [15-17]. Therefore, it is also 
highly desirable to find approaches to determine Uc that 
could be more easily implemented in experiments. 

In this work, we address the two issues mentioned 
above, namely, how to accurately determine Uc within 
computational approaches and in experimental studies. 
First, we apply a recently proposed scaling approach for 
the gap [18] to obtain the critical point in the BHM with 
high accuracy and at fillings n = 1, 2, 3. Second, we in- 
vestigate the behavior of the fidelity susceptibility across 
these transitions. The fidelity susceptibility, a quantity 
motivated from the field of quantum information, has re- 
cently attracted much attention as a means of identifying 
the presence of quantum phase transitions even if the na- 
ture of the involved phases is not known [19-28]. Finally, 
we discuss how to determine Uc by studying the center of 
mass motion during Bloch oscillations, which occur after 
tilting the lattice [29, 30]. This is something that can be 
easily implemented in ultracold gases experiments. 

All calculations are done utilizing the DMRG and, out- 
of-equilibrium, the Krylov- variant of the adaptive time- 
dependent DMRG (t-DMRG) [8]. For ground state cal- 
culations, we perform 10 sweeps and keep up to m = 1000 
density matrix eigenstates. In order to ensure the high 
accuracy needed for the considerations below, we trun- 
cate the local Hilbert space at n -I- 5 bosons (where 
n = 1, 2 or 3). The ground state energies obtained are 
converged in most cases with an absolute accuracy of 
10~^ or better. The computation of the fidelity suscepti- 
bility is, however, more demanding, and we are restricted 
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to smaller systems. In order to reach the necessary accu- 
racy in the overlap of the two wave functions involved, we 
use m < 4000. To study the out-of-equilibrium dynamics 
under a tilt, we truncate the local Hilbert space at n + 4 
bosons and keep up to to = 2500 states in the course of 
the time evolution while using a time step of At = 0.01. 

Scaling analysis of the gap — The phase transition from 
a Mott-insulator to a superfluid in the ID BHM at com- 
mensurate fillings is known to be of BKT type [31]. 
Hence, it is accompanied by the exponential closing of 
the single-particle gap Eg ^ exp{—b/y/U — Uc) {b is a pa- 
rameter which is independent of U). As a consequence of 
its exponential behavior, a direct study of the transition 
by computing the single-particle gap for finite systems is 
plagued by finite-size effects. This problem can be over- 
come by a scaling analysis of the gap, for which we follow 
the approach in Ref. [18], briefly described below. 

The method is based on the following ansatz for the 
scaling of the gap in the vicinity of the phase transition, 
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where is a scaling function, C is an unknown con- 
stant to be determined, and L is the system size. We 
emphasize two aspects of this scaling ansatz: First, it 
contains the logarithmic corrections that are typical for 
Eg (L) at the BKT transition [32, 33]. Second, it resem- 
bles the relation for the resistance (which also vanishes 
exponentially) in the charge-unbinding transition of the 
two-dimensional classical Coulomb gas, which is also of 
BKT type [34]. At the critical point, and in its vicin- 
ity within the superfluid region, one expects the values 
of F{^/L) to be system-size independent because of the 
divergence of the correlation length. Hence, the data for 
the rescaled gap E*g (L) = LEg (L) [l + l/{2\nL + C)] 
for different system sizes L will be independent of L in 
this region. Furthermore, the curves E* (L) vs ^/L for 
several values of L and U should collapse onto a unique 
curve representing F. Equivalently, one can reformulate 
the relation in Eq. (2) by taking the logarithm of the ar- 
gument of F and considering a different function / with 
argument = InL — In^. 

We determine the critical point by adjusting the pa- 
rameters Uc, b, and C. In the procedure, we look for the 
best collapse of the curves E* (L) vs a;^ for different val- 
ues of U and L. This is done by representing the function 
/ with a selected high-degree polynomial (eighth degree 
in our case) such that the results are independent of the 
degree. Such polynomial is flt on a dense grid of values 
of Uc, b, and C, to the calculated values of E* (L) and 
XL. The quality of the fit is assessed by computing the 
sum of squared residuals (S), which defines the function 
S{Uc, b, C). Uc is then obtained from the set of parame- 
ters Uc, b, and C which minimizes S{Uc,b,C). 

We have applied this procedure to integer filled chains 
with n = 1, 2 and 3. We find that, in these three cases, 
the minimum of S{Uc, b, C) is obtained for arbitrarily 
large values of C . This means that logarithmic correc- 
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FIG. 1. (Color online) (a) Contour plot of the sum of squared 
residuals S{b, Uc) for n = 1. The white arrow signals the 
location of the minimum value of S. The white lines are 
equally spaced contour lines where S is constant, (b) Best 
collapse of the data for Eg (L) vs xl corresponding to Uc = 
3.279, b — 5.2, and C — >■ oo. The inset shows the rescaled gap 
vs. U. A similar analysis for n = 2 [n = 3] is presented in 
panels (c) and (d) [(e) and (f)]. 



tions to the scaling of the gap, in the form (2), do not 
play a role in the determination of the critical point. This 
is to be contrasted with the t-V-V model in Ref. [18], 
where C was found to be finite in all transitions analyzed. 
In Fig. 1(a), we present a density plot corresponding to 
S{Uc,b, oo) for n = 1, which exhibits a clear minimum at 
Uc = 3.279 ± 0.001, 6 = 5.2 ± 0.1. The error bars are es- 
timated by repeating the minimization procedure adding 
and subtracting to the gap the error of the energy (as- 
sumed to be = 10~^). Furthermore, the sensitivity 
of the results to the selection of the interval of values of U 
used in the fit is also included in the error bars such that 
our results are independent of that choice. Correspond- 
ing to the set of parameters that minimizes S{U, b, C), in 
Fig. 1(b), we plot E* (L) vs xl- The data is clearly seen 
to collapse to a single curve representing the function /. 
In the inset, the curves for the rescaled gap correspond- 
ing to different system sizes are seen to merge around 
the critical value Uc- Our result for the critical point is 
in agreement with recent calculations by other authors 
(see, e.g., Ref. [11] where Uc ~ 3.27). Analogously, we 
find Uc = 5.587 ± 0.001 for n = 2 [Figs. 1(c), and 1(d)], 
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FIG. 2. (Color online) Fidelity susceptibility for different sys- 
tem sizes at integer filling n = 1 (red), n = 2 (green) and 
n = 3 (blue). The plot shows data for L — 40 (square), 
1/ = 80 (circle) and L = 120 (triangle), the solid lines are 
spline interpolations and serve as a guide to the eye. The 
black lines (diamonds) are the result of a finite-size extrapo- 
lation using a quadratic fit. The vertical dotted lines indicate 
the position of the quantum critical points obtained using the 
scaling analysis of the gap described in the text. 



and Uc = 7.876 ± 0.002 forn = 3 [Figs. 1(e), and 1(f)]. 
We use these results in the following to benchmark alter- 
native approaches for estimating the transition point. 

Fidelity Susceptibility — The fidelity susceptibility (FS) 
X for the ground state of the system \iPq) is defined as 



xiu) 
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and is also known as the fidelity metric. For generic 
second order phase transitions, x is expected to diverge 
in the thermodynamic limit (TL) [19, 21-24], and it has 
been found to exhibit clear signatures of such transitions 
already for rather small system sizes, where a maximum 
of X is seen near the transition point [25-28]. 

In Fig. 2, we show the fidelity susceptibility for the 
BHM at fillings n = 1, 2, 3, for system with L = 
40, 80, 120, and for on-site interactions up to [/ = 8. 
For n = 1 and n = 2, x exhibits clear maxima for val- 
ues of U greater than Uc computed from the scaling of 
the gap. For n = 3, the maximum is expected to be be- 
yond the values of U studied here. Hence, indications for 
the existence of a phase transition are obtained already 
for small systems. Interestingly, and also of relevance to 
the x's calculated here, recent works have proposed that 
a minimum of the FS may signal the quantum critical 
point [35, 36]. This was argued to be possible because, 
depending on the scaling dimensions of the system, the 
FS can be finite at a critical point [21]. In Fig. 2, one can 
indeed see that a minimum of x also occurs close to the 
critical point. Furthermore, between the maximum and 
the minimum, there is a point at which all values of x 
seem to be independent of L for the system sizes treated. 
As seen in Fig. 2, for n = 1 and n — 2, this occurs at 
values of U greater than Uc computed from the scaling of 
the gap. We applied different extrapolation schemes for 



X and did not obtain results consistent with those from 
the scaling of the gap: For example, in Fig. 2 we show 
the outcome of the simplest approach in which the ex- 
trapolation to the thermodynamic limit is attained using 
a second order polynomial. Neither the position of the 
maximum, nor the one of the minimum or of the cross- 
ing point are in agreement with the values of Uc obtained 
from the scaling of the gap. As pointed out in previous 
studies (see, e.g., Ref. [24] for an analysis of the ID Fermi 
Hubbard model) , the divergence of x could be extremely 
slow and very large system sizes (as well as a more elabo- 
rate finite-size-scaling ansatz) may be required to resolve 
the critical point. Further studies are needed to address 
this problem in the BHM. 

Center of mass motion after tilting the lattice — In or- 
der to determine the critical point in experiments with 
ultracold bosons in optical lattices, we propose to follow 
a recent proposal that uses Bloch oscillations [29]. The 
idea is to apply an external field 
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at time i > 0, to a system that is initially in its ground 
state {n = ior t < 0). Such a set up can be realized 
in optical lattice experiments by, e.g., tilting the lattice. 
One can then study the center of mass motion (COM) 



xcuit) = ^J2j{nj){t), 
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where TV is the total number of particles, at times t>Q. 
In previous studies, in a variant of the t-J model at low 
filling [29] and in an effective Ising model in a transverse 
field [30], it was observed that the amplitude of the COM 
exhibited signatures of the quantum phase transition at 
the critical point. In the t-J like model [29], because 
of the formation of pairs, the transition from a metallic 
to a gapped superconducting phase was visible by both 
a rapid decrease of the amplitude and by a doubling of 
the frequency of the Bloch oscillations. In the Ising like 
model [30] , the amplitude was found to be maximal at the 
transition point. From the experimental point of view, 
Bloch oscillations have been, e.g., used to investigate 
Dirac points on hexagonal optical lattices [37], as well 
as to study low-frequency breathing modes in elongated 
Fermi gases [38]. Here, we investigate what happens in 
the BHM and, at the same time, analyze the simpler (in- 
tegrable) case of spinless fermions with nearest-neighbor 
interaction V , 



H.c 



(6) 



This model is exactly solvable and, at half filling, also 
possesses a BKT transition from a Luttinger liquid (LL) 
to a charge density wave (CDW) insulator at V/ J — 2 
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FIG. 4. (Color online) Amplitude of the first oscillation of 
the COM after tilting the lattice, (a) Spinless fermions at 
half fiUing (L = 20 and L = 40) and ^ = 1. (b) BHM at 
filling n = 1 (L = 14 and L = 20) and H = 1. The inset 
shows results for L = 14 at filling n = 2. 



FIG. 3. (Color online) COM for (a) spinless fermions {L = 20, 
f2 = 1) after tilting the lattice for the indicated values oiV = 
0, . . . , 3. (b) the BHM {L ^14,0. = 1) after tihing the lattice 
for the indicated values of U — 1, ... ,5. The red stars and 
blue circles denote maxima and minima of the oscillations, 
respectively. 



[2, 18]. In what follows, we set J = 1 and h = 1, so that 
V is given in units of J and time t in units of h/ J. 

In Fig. 3(a), we display the COM of a half-filled chain 
of spinless fermions with L — 20 and different values 
of V on both sides of the LL to CDW transition. As 
seen there, the amplitude of the oscillations decreases 
and the damping rate increases when increasing V and, 
on the CDW side, no oscillations can be resolved. This 
is reminiscent of the behavior of a harmonic oscillator 
which moves freely {V = 0), damped {0 < V < 2) and 
overdamped ^ 2). A qualitatively similar behavior is 
observed in the COM of the BHM at filling n = 1 and 
for 1 < < 5 [Fig. 3(b)]. There, finite values of U < 3.5 
lead to damped oscillations, and only the first oscillation 
can be resolved on the time scale of our simulations. For 
U > 3.5, overdamped behavior sets in and no oscillations 
can be identified. 

To gain a better understanding of the evolution of the 
Bloch oscillations as interactions are increased, in Fig. 4, 
we display the amplitude (defined as the difference 
between the first maximum and the first minimum), 
for spinless fermions vs V [Fig. 4(a)] and for the BHM 
vs U [Fig. 4(b)]. For spinless fermions, it can be seen 
that the amplitude of the oscillations at the critical 
point and above {V > 2) is very small and decreases 
with increasing system size. The results for the BHM 
are qualitatively similar. The region of U at which the 
amplitude of the Bloch oscillations is seen to vanish, 
3.0 < U < 3.5, contains the value obtained from the 
scaling analysis of the gap Uc ~ 3.279. Therefore, for 
both systems, the BKT transition leads to comparable 
behavior, and the study of Bloch oscillations provides a 
good experimental estimate of the critical value for the 



We also studied the COM for n = 2 and n = 3. Since 
computations become increasingly demanding with in- 
creasing filling, only smaller lattice sizes could be stud- 
ied in those cases. In addition, the numerical values of 
the amplitude become significantly smaller. We therefore 
find that theoretical estimates for Uc from the behavior 
of the Bloch oscillations become less accurate with in- 
creasing n. In the inset in Fig. 4(b), we show results for 
ri = 2 and L — 14, where one can see that the ampli- 
tude of the Bloch oscillations vanishes for 4.5 < U < 5, 
while the theoretical prediction is Uc = 5.59. It would be 
interesting to experimentally study Bloch oscillations in 
larger lattice sizes and see if the worsening of the predic- 
tion is due to finite size effects or due to a worsening of 
this approach with increasing filling. 

Conclusions and Outlook — We followed three ap- 
proaches to study quantum critical behavior in the one- 
dimensional BHM at integer fillings. By means of a scal- 
ing analysis of the gap, we obtained accurate values of 
Uc for the superfluid to Mott insulator transition at fill- 
ings n — 1, 2, 3. The fidelity susceptibility was shown 
to exhibit signatures of the phase transitions for finite 
systems, but the results for this quantity did not allow 
us to improve on the values of Uc obtained from the scal- 
ing of the gap. Finally, we showed that the study of 
Bloch oscillations in experiments can help locating the 
critical values for the superfluid to Mott insulator transi- 
tion. The latter approach could potentially be used also 
in experiments in higher dimensions. 
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